Anomalous spin and charge dynamics of the t-J model at low doping 
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We present an exact diagonalization study of the dynamical spin and density correlation functions 
in small clusters of 2D t~J model for hole dopings <25%. Both correlation functions show a 
remarkably regular, but very different scaling with both hole density ph and parameters t and J: 
the density correlation function is most consistent with that of condensed Bosons in a band of width 
the spin correlation function with that of fermions in a band of width ~ J. We show that the 
familiar spin bag scenario explains these results in a natural way. 
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The identification of a simple 'effective theory' for the 
doped Mott-Hubbard insulator which is capable of re- 
solving the numerous anomalies of cuprate superconduc- 
tors remains an intriguing problem. The simplest model 
which may be expected to contain the key features of 
these systems is the t—J model: 
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The Si are the electronic spin operators, c\ ^ — c\ — 
Tii^^a-) and the sum over < i,j > stands for a summation 
over all pairs of nearest neighbors. 

For hole densities ph>0.3 the ground state of the model 
seems to represent a fairly conventional Fermi liquid with 
a particle-hole like spin and density excitation spectrum 
10,^. In this work we show that the situation is drasti- 
cally different for p/i<0.25: here the density excitations 
roughly resemble those of condensed Bosons with charac- 
teristic energy t, the spin excitations still are consistent 
with Fermions, the characteristic energy, however, now 
is J. We show that these results are naturally explained 
within the familiar spin bag or string scenario. 
For the standard 16 and 18-site clusters we used the Lanc- 
zos algorithm to compute the dynamical spin (SCF) and 
density (DCF) correlation functions: 
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Here j'l'o) (Eq) denotes the ground state wave func- 
tion (ground state energy), for the operator Oa{q) we 
choose the Fourier transform of either density operator 
(ni_|-|-ni^l) (a=d) or spin operator (jT-i,]-— fii.j,) {a=s). 
Fig. shows the DCF divided by the number of holes, rih, 
for n/i=l, 2, 3, 4. Except for the w — > parts at (7r/2, 0) 
and (7r/3,7r/3) this is a universal function. Deviations 
are strongest for a single hole, but even there the char- 
acteristic shape is already present. Fig. thus first of all 
demonstrates a high degree of continuity over the entire 
range of dopings considered; quite obviously the essential 



physics is already realized for a single hole in an antifer- 
romagnet. Next, the scaling of the spectra with hole 
concentration over their entire width is clearly inconsis- 
tent with particle-hole excitations in a Fermion system: 
at least the high-energy parts, where transitions from 
deep below to high above the Fermi energy Ep would 
contribute, should be unaS'ected by a change of particle 
density. Instead, the DCF could be modelled roughly by 
Bosons of density ph which are condensed into the lowest 
state of the noninteracting band with the free electron 
dispersion Cg; such a system would have a DCF of the 
form Cd{q,i^)—ph ■ S{uj — eq). That the characteristic 
energy of the DCF indeed is t is seen in Fig. , which 
compares the DCF for different t/J and demonstrates 
that the positions of the dominant 'peaks' remain unaf- 
fected by a change of J. Fig. also shows that the weight 
of the LU ^ part decreases with t/J, whereas the weight 
of the remainder other parts increases; together with the 
deviations in the scaling behaviour with rih this suggests 
a two-component interpretation of the DCF. 
We proceed to the SCF, shown in Fig. for n;i=2,3,4. 
The spectra for 3 holes thereby are averaged over the 4 
degenerate ground states (the 3 hole ground states have 
momentum (27r/3, 0) in the 18 and {tt/2,tt/2) in the 16- 
site clusters). With the exception of (7r,7r), the spectra 
are fairly independent of doping, changes occur mostly 
at the low-energy end of some of the spectra. Such 
a dependence on particle number is reminiscent of the 
particle-hole excitations in a Fermion system: the spec- 
tral weight of a given transition is expressed in terms of 
the Fermionic occupation numbers rik as nfc+q(l — rik) 
and is affected only if either initial or final state cross 
the chemical potential. Here it should be noted that due 
to the coarse fc-mesh of the 2D clusters (as well as the 
averaging process for 3 holes) one can not expect a well 
defined kp which scales continuously with particle den- 
sity. Instead a countinuous increase of the hole occupa- 
tion numbers at the 'Fermi momenta' and hence a 
continuous increase/decrease of peak intensities will be 
realized. As for the energy scale Fig. shows the SCF for 
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various t/ J. When frequencies are measured in units of 
J, the peak positions obviously remain largely unaffected 
by a change of i, so that the characteristic energy of the 
SCF is J. 

Summarizing the numerical results one may say that both 
SCF and DCF scale with nu and t/J in remarkably reg- 
ular but completely different ways. Let us stress for clar- 
ity that the different doping dependence of the integrated 
weight with nonvanishing momentum transfer is trivial. 
By the /-sum rule the integrated spectral weight of both 
correlation functions equals the number of electrons, rig. 
For the DCF this value is almost exhausted by the q=Q 
spectrum, which contributes ne{ne/N) (with N the num- 
ber of sites). For the SCF this contribution either van- 
ishes or is a tiny 1/iV, so that the integrated spectral 
weight for finite q is ~rte for the SCF but for the 
DCF. This argument, however does not explain e.g. the 
scaling with Uh of the entire DCF spectra. This is not 
only inconsistent with the particle hole excitations in a 
Fermi liquid but also with the Luttinger liquid realized 
in the ID model 0J^. The low doping phase in 2D thus 
obviously exhibits qualitatively new physics. 
To gain additional insight. Fig. compares the SCF for 
t/J— 2.5 and t=0, i.e., mobile and static holes. For static 
holes there is a band of diffuse excitations, which roughly 
follow the characteristic spin wave dispersion; their dif- 
fuse nature is obviously due to the scattering from the 
holes, which in this case merely act as impurities. For 
mobile holes spin-wave excitations near (tt, tt) seem to 
persist with reduced spectral weight and a wider gap at 
(tt, tt); the qualitatively new feature are low energy peaks 
throughout the Brillouin zone which may naturally be 
associated with particle-hole transitions in the coherent 
band for mobile holes. This suggests a two-component in- 
terpretation for the SCF, namely low-energy particle-hole 
excitations plus spin- wave like excitations near (7r,7r). 
This is also supported by the opposite t/ J-dependence 
of these two components (Fig. ) : whereas the weight of 
the 'spin wave' parts decreases with increasing t/J (con- 
sistent with the more efficient degradation of antifcrro- 
magnetic correlations by more mobile holes) that of the 
particle-hole-excitations increases or remains unchanged. 
Further information is obtained by studying the impact 
of Os and Od on the electronic momentum distribution 
(EMD) n(fe) = {Cka^k,a)- To that end we compute 
first the ground state EMD, no(fe), and second the EMD 
nq{k) for the state {l/n)Oa{q)\GS) (where n is chosen 
to normalize the state to 1). Table | shows the dif- 
ference Aq{k)—nq{k)—no{k) for q — (tt, tt), as well as 
nshiftiq) = J2k \^qif^)\- The latter quantity may be in- 
terpreted as the number of electrons shifted in fc-space, 
its value for free particles is 1. Table | reveals a clear 
difference between density and spin operator: whereas 
Od induces a substantial shift of electrons in k space and 
always has n shift close to its free-particle value of 1, Os 
affects n{k) to a much lesser degree. Via the kinetic en- 



ergy sum rule 

-Bfc,„ = 2^efcn(fc), (1) 
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(with efe = —2t{cos{kx) + cos(fcy))), this result on one 
hand immediately explains the different energy scales 
of SCF and DCF: Od substantially changes Ekm'^t, Os 
leaves it essentially unchanged. On the other hand, this 
only complicates the puzzle about the doping dependence 
of the DCF: quite obviously the electrons react to the 
density operator as if they were free particles, hence one 
naturally would expect a Fermion-like particle-hole spec- 
trum for the DCF, in contrast to the numerical result. 
We now want to show that all results obtained so far can 
be resolved in a simple and natural way if one adopts the 
familiar string or spin bag picture. Thereby we consider 
holes moving in (and hence coupled to) a 'background' 
of antiferromagnetically correlated spins. Our key as- 
sumption is that the spin background carries excitations 
which are independent of the hole system, so that their 
momentum is not 'visible' in the EMD (compare Table |) 
and hence does not increase the kinetic energy. Natural 
candiates are the remnants of the short-wavelength spin 
waves, as would be suggested by Fig.. We next assume 
that the relevant hole states are described by operators 
of the type 

Cu,k,^ = a,(fc)cfe,| + ^ ^ (fc, k')ck',aSCk' 

+ E E ir'^"'\k,k',qycy.,,,s'(sCk'-q 

k' ,q (T,(T' ,fT" 

+ ..., (2) 

where denotes the electronic spin operator. Equa- 
tion (||) describes a hole which has transferred a part 
of its momentum to a variable number of spin excita- 
tions; it has many degrees of freedom (as reflected by 
the many parameters a, /3, 7 . . .) so that there will be a 
large number of bands, which presumably form the ex- 
tended incoherent continua present in the single particle 
spectral function The bands are labelled by the in- 
dex v, denotes the 'quasiparticle band' |^ split off 
from the bottom of the continuum. For this quasiparticle 
band the validity of the spin bag description has been ver- 
ified previously by explicit numerical check The 
driving force behind the momentum transfer to spin ex- 
citations is gain in kinetic energy: scattering a hole e.g. 
from fe=(7r/2, 7r/2) to fc'=(7r,7r) and transferring the ex- 
cess momentum to a spin excitation reduces the kinetic 
energy ^t, but requires only exchange energy ~ J<i. We 
can conclude that in the low lying spin bag states the 
bare hole predominantly occupies momenta near (7r,7r): 
via (|l|) the depletion of n{k) on these momenta lowers 
the kinetic energy most efficiently. This is corroborated 
by the numerical result |jll[ that the addition of e.g. a 
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single hole with total momentum (7r/2, 7r/2) reduces n(fc) 
strongest near (tt, tt) . A rough estimate for the degree of 
admixture of spin excitations is the quasiparticle weight 
Z: the simplest estimate would be ^=|ao(fcF)|^; -^—0.3 
[ p2[ thus suggest a rather strong admixture of spin fluc- 
tuations. 

We now assume that in the Uh hole-ground state the 
holes occupy the Uh lowest states of the ^=0 band, i.e. 
rigid-band filling of the quasiparticle band; this is con- 
sistent with numerical results Since the dispersion 
of the j/=0 band strictly scales with J ||] , we may con- 
clude from (|l|) that the distribution of the bare holes in 
fc space is essentially the same for all v=Q states. In 
other words, the coefficients /?, 7 in (||) depend strongly 
on the 'bare hole' momentum fc' but only weakly on the 
total momentum fc. We can conclude that when holes 
are filled into the v—Q band, the hole occupation of e.g. 
(tt, tt) should be times that for a single hole; that this 
indeed is the case is shown in Table O which gives the 

hole occupation, ''^(fc) = (l/2) •^^^("^fc-o''^ o-)' ^"-"^ momenta 
near (7r,7r) in the ground state with different Uh. 
With this picture of the ground state in mind, we can dis- 
tinguish two types of excitations: there can be 'particle- 
hole excitations' within the v — Q band or an excitation 
of the 'internal degrees of freedom' of a spin bag, where 
the final state has v' An important point is that Os 
and Od transfer momentum to a spin bag in a very differ- 
ent way: Os changes the momentum of the bag by adding 
(or removing) a spin excitation, whereas Oa can transfer 
its momentum only to the bare hole itself. Os conse- 
quently does not change the distribution of holes in fc- 
space appreciably whereas Od necessarily must do so (see 
Table |). Hence, Os is essentially limited to transitions 
within the i^=0 band, so that its spectrum is Fermionic 
and scales with the quasiparticle bandwidth J. On the 
other hand, Od induces transitions from the v=Q band 
to 'bands' in the continuum, with an excitation energy 
~t which usually far exceeds the bandwidth ^J. Then, 
particle-hole transitions of Fermions between bands of 
width W ^ which are separated in energy by E~^W ^ will 
give the contribution C(lo) ~ Pf5{L0 — E), to the cor- 
relation function, where Pf is the density of Fermions 
in the lower (partially occupied) band. This immedi- 
ately explains the apparently Bosonic doping dependence 
of the DCF, which thus originates from the possibility 
to excite high-energy degrees of freedom of the struc- 
tured spin bag quasiparticles. Only for small momentum 
transfer the increase in kinetic energy is small, so that 
Od also can generate particle-hole transitons in the i/=0 
band, where restrictions due to the Pauli principle apply, 
hence the deviations from the scaling of the DCF with 
n/i in the low energy region. The spin bag scenario thus 
provides a natural explanation for the unusual scaling 
behaviour of the different correlation functions; it is sup- 
ported by detailed consistency with a substantial body 
of numerical evidence. For completeness we note that 



slave-boson mean-field theories describe the ground 
state of the t—J model as a product of condensed Bosons 
in a band of width and Fermions in a band of width 
r-^AJ. With the additional assumption that the density 
operator acts only on the bosons, the spin operator only 
on the Fermions, this ground state clearly would have an 
excitation spectrum which is consistent with the numeri- 
cal results; the justification of this assumption as well as 
the possible agreement with details of the cluster results 
remains to be clarified. 

In summary, we have studied the dynamical spin and 
density correlation function for the 2D t—J model near 
half-filing. Whereas these correlation functions should be 
closely related in a Fermi liquid, we found them to differ 
substantially for this strong correlation model: the den- 
sity correlation function has a Boson-like dependence on 
the hole density and the hopping integral t as its charac- 
teristic energy scale, the doping dependence of the spin 
correlation function is consistent with Fermions and it 
has the exchange constant J as energy scale. While the 
remarkably systematic scaling of the correlation func- 
tions suggests the existence of a simple 'effective theory' 
for the excitation spectrum, the familiar particle-hole pic- 
ture thus is clearly insufficient. The familiar spin bag 
picture then provides a promising framework for such an 
effective theory: the strong dressing of the hole with spin 
fluctuations and the resulting complex internal structure 
of the spin bag-type quasiparticles lead to a qualitatively 
new type of excitations, namely the excitation of inter- 
nal degrees of freedom of the quasiparticles. Spin and 
density operator differ markedly in their ability to ex- 
cite the various degrees of freedom of the spin bag liquid, 
hence their very different spectra. While details need to 
be worked out, it seems obvious that the existence of 
such new types of excitations, as well as the apparently 
very different response of the spin bag liquid to 'spin- 
like' and 'charge-like' perturbations may lead to experi- 
mentally observable anomalies e.g. in high-temperature 
superconductors. 

Financial support of R. E. by the Japan Society for the 
Promotion of Science is most gratefully acknowledged. 

DCF divided by rih for various rih [t/J = 2.5, 
Lorentzian broadening e=0.2i). 

DCF for 2 holes in 16-sites for different t/J (e=0.2t). 

SCF for different rih] the (7r,7r) spectra are multiplied 
by 0.2 (VJ = 2.5, e=0.5J). 

SCF for 2 holes in 16 sites for different t/J; the (tt, tt) 
spectra are multplied by 0.2, e=0.2J. 

SCF for 2 mobile and static holes. The (tt, tt) spectra 
are multipfied by 0.2, e=0.5J. 
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TABLE I. (a) A 
q={TT,n) (16 sites, 2 holes, t/J 
same system. 



q(fc) for spin and charge operator with 



2.5). (b) Ushiftiq) for the 
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TABLE IL Hole occupation numbers n(fe) divided by 
for different k and Uh (16 sites, t/J = 2.5). 
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Figure 2 
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Figure 3 
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Figure 4 
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Figure 5 
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